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Abstract 

We study theoretically the surface of time-reversal-symmetric, noncentrosymmetric supercon- 
ductor with mixed singlet and triplet order parameters. A pair of counterpropagating subgap 
quasiparticle surface bound states with opposite spin projections are obtained in the nontrivial Z2 
case where the triplet component is larger than the singlet one, contributing to a spin current with 
out-of-plane spin projection. In contrast to the pure p-wave cases, these subgap states do not have 
fixed spin projections, which however depend on the momenta along the surface. In the trivial Z2 
case where the singlet order parameter is larger, no subgap surface bound states show up. In both 
cases, there is also a finite contribution to the spin current from the continuum states with energies 
between the two gaps. The method for obtaining the quasiclassical Green's functions associated 
with the noncentrosymmetric superconductors is also presented. 
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I. INTRODUCTION 



How to manipulate spin in condensed-matter system has been the main challenge for 
both experimentalists and theorists in this community. Recently, the predictiona^"- and 
observations^"- of topological insulators (TI) with time- reversal symmetry (TRS), such as 
HgTe/ (Hg,Ce)Te semiconductor wells,- Bii_a,Sba:,- Bi2Se3,^ and Bi2Te3,- inspire a great deal 
of interest in both the application and the fundamental research ends . The main feature of 
a TI is a pair of counter-propagating edge states with opposite spin projections developed 
out of a gapped band structure, very much like the edge states in the integer Quantum 
Hall case except that the TRS is broken in the latter case. In the language of homotopy, 
the mapping from the momentum space to the Hamiltonian can be smoothly deformed into 
either one of the two distinct elements in the so-called Z2 class associated with the trivial 
and nontrivial TI's. This pair of surface bound states contribute to a spin current near the 
surface of a TL- 

The concept of TI can also be generalized to superconductors (SC). A simple example 
is the (2D) p-wave superconductor with its order parameter given by ci = kxV — kyX. We 
can see explicitly from the order parameter oc [{k^ — iky)\ tt) + {k^ + iky)\ H)] that the 
Cooper pairs with down(up) spins have a counterclockwise(clockwise) motion. This state 
is clearly time- reversal symmetric. Each species is in an axial state, considered in, e.g., 
Ref.-*^, with the phase of the Cooper pair wavefunction advances (decreases) by 2it when 
the angle of the momentum direction advances by the same angle. Considering the surface 
of this superconductor adjacent to vacuum, the incident and reflected quasiparticles see an 
order parameter with a different phase factor, analogous to the case of a Josephson junction. 
A pair of surface bound states with opposite spin projections propagate in opposite direc- 
tions, generating no charge but finite spin current, as in the case of a topological insulator. 
Moreover, straightforward generalization of the results of Ref.-^^ implies that a singly quan- 
tized vortex of this superconductor possesses a pair of zero-energy Majorana states inside 
its vortex core. In contrast, the order parameter of an s-wave singlet superconductor do not 
change sign over the Fermi surface. No surface bound states are topologically required and 
the superconductor thus belongs to trivial class. 

For noncentrosymmetric superconductors (NcSC), such as the compounds CePtaSi,— 
Li2Pt3B,^^ and CeRhSis,- and the 2D electron gas between two insulating layers,— complica- 
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tion arises due to the presence of parity-broken spin-orbital interaction and the singlet-triplet 
mixed superconducting order parameters.— For example, for a system with no up-down re- 
flection symmetry, the s-wave singlet order parameter and the triplet order parameters 
Ap with the d-vector given by (fc^^y — kyx) would naturally mix since they are of the same 
symmetry. Alternatively, the two Fermi surfaces with opposite helicities can be associated 
with two different superconducting g result of the broken inversion symmetry.— 

Nevertheless, the topological classifications still can be deduced from the existence of zero- 
energy vortex bound states^ or the surface bound states.— These studies reveal that only 
the relative signs of the pairing terms on the opposite-helicity bands matter for the topology, 
which can also be shown by more explicit topological arguments.— i^i NcSC with opposite 
signs of pairing on the two bands resembles the pure p-wave triplet superconductor in that a 
pair of topologically protected zero energy or surface bound states reside within the vortex 
core or at the surface, respectively, while these topological bound states do not exist at all 
in the NcSC with same sign of pairing, which resembles the pure s-wave singlet SC. 

In this paper we shall consider in more detail the surface of a NcSC, as shown in FigH], as 
a function of the singlet A^ and triplet Ap order parameters. For simplicity, here we would 
not consider the dispersion along the z axis (or effectively 2D) nor the splitting of the fermi 
surfaces due to spin orbit interactions. The absence of spin-orbital coupling may in fact 
correspond to a narrowly accessible physical regime in mixed-parity superconductivity.— 
Our simplification however does not affect the main physics that we would like to explain 
in this paper. Comments on this aspect will be given later. 

For the present geometry, there is no xy plane reflection symmetry, whereas the time- 
reversal symmetry and the xz plane reflection symmetry are left intact. Following the same 
symmetry arguments as in our previous paper,— the magnetization along any directions 
automatically vanish due to the time- reversal invariance. However, the spin currents Jy, Jy, 
and J| can in principle be allowed. We shall examine whether these spin currents are finite. 

In the pure p-wave case (|As|=0), the bound states have their spins quantized parallel or 
antiparallel to the z-axes. Only Jy is non-zero.—'^i We would examine what happen to these 
states when A^ is finite, especially their spin directions. We shall see that the spins are no 
longer polarized along z in the general case. In the context of topological superconductor, 
the parity- mixed order parameters with |Ap| > |As| belong to the nontrivial Z2 class,—""— 
and the pair of surface bound states with opposite spin projections propagating in opposite 



3 



directions still exist. These surface bound states can in principle be detected by tunneling 
conductance, which has now been investigated in great detail theoretically.— ^^^i^ We shall 
however concentrate on the spin currents in this paper. These surface states can still generate 
the spin currents Jy (and in principle also J^', but see below). However, as we shall see, 
this is not the only contribution. For the case of |As| > |Ap|, the superconductor belongs 
to the trivial category and it is expected that the surface bound states do not survive. An 
immediate question is, at the transition point for |Ap| = |As|, whether the spin current 
drops to zero abruptly, signaling the transition of topology, or it is smoothly decreasing 
towards zero, coinciding with the pure s-wave case only when Ap vanishes. In the following, 
we shall consider the above questions by evaluating quantities such as the momentum- and 
spin- resolved density of states and the spin current of the NcSC with varying ratio |Ap|/|As|. 
We find that nonzero spin current Jy is dictated only by the broken symmetry and can be 
finite for both the topological trivial and non-trivial superconductors. The contribution to 
this spin current turns out to arise also from non-topological continuum states with energies 
between the two gaps |Ap ± As|. These states are present in both the topological trivial 
and non-trivial superconductors, and they are not required by topology. (We recall here the 
analogous situation that a finite spin Hall conductivity is possible from the Kubo formula:^ 
yet the model can belong to the topologically trivial Z2 class^). Lastly, though in principle 
Jy, and jy are allowed by symmetry, we found that they vanish within our calculations. 

II. QUASI-CLASSICAL GREEN'S FUNCTION ASSOCIATED WITH NCSC 

Now we use the quasiclassical Green's functions to investigate the surface of a clean NcSC. 
Here we shall employ the so-called exploding and decaying tricks,—"^ which is related to the 
projector formalism initiated by Shelankov.'^^ This method is different from the approach 
which employs the Riccati transformation.-S^ The Matsubara Green's function g(k, e„, R) in 
spin and particle-hole space satisfies the Eilenberger equation,— 



[ze„r3 - A, g] ivf ■ Vng = , 
with the normalization condition 



(1) 



Here e„ and k denote the Matsubara frequency and momentum direction associated with 
the quasiparticles, respectively. R represents the spatial position. The set of matrix {1, r} is 
used in the particle-hole sector while {1, a} serve for the usual spins. In this representation, 
the 4x4 pairing order parameter can be written as 

where the 2x2 matrix A=[As + Apd(k) ■ (y]{iay) in the usual spin representation. Here • 
and • are to denote the 4x4 and 2x2 matrices, respectively. The triplet order parameter 
we consider is of Rashba form d = {—ky, k^). As and Ap can be both taken real since we 
assumed TRS. Without loss of generality, they will both be assumed positive. There are two 
energy gaps, | Aj, ± | , associated with the quasiparticles in the superconducting states.— 
One should also note that we have not included the spin-orbital coupling term which may 
arise from the lack of inversion symmetry in the normal state. Therefore, only one Fermi 
surface and its associated Fermi velocity vj are needed. Qualitative effects of including the 
spin-orbital coupling will be discussed later in the paper. For simplicity, we shall also ignore 
the spatial variations of the order parameters As^p. 

We consider the geometry shown in Fig. [1] The incoming and reflecting quasiparticles 
have the momentum k and k, respectively. Here k^ > and < 0. We label positions along 
the quasiparticle path consisting of each pairs of k and k by u, with u < {u > 0) labels 
the part for k (k). The quasiclassical Green's function at the surface g(0) can be obtained 
in terms of the decaying and exploding solutions to eq. ([I]) (c.f. Ref.— ) 

g(0) = -m[i(k),B(k)]{A(k),B(k)}-i , (4) 

where the matrix A = a'^^ + a and B = 6++ -|- b . Below the 4 matrix solutions a"^"^ 
and are explicitly shown, 

a^^ = (^^^^)H|A±|V3 - z(a± + e„)(A, + Apd ■ a)ayr+ (5) 

+i{a± - en)ay{Al + A*/ ■ a)r_] , 

b^^ = (^^^^)HA±|V3 - t{a± - e„)(A, + Apd-a)ayT+ (6) 

+i{a± + en)(Jy{Al + A*pd ■ a)r_] , 
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where the notation we used for the three 4x4 matrix a = {crx,(TyT3,az) should be noticed. 



The parameters a± = ye^ + |A±p, where A± = Ag ± Ap. More details can be found in 
Appendix El 

III. SURFACE BOUND STATES AND SPIN CURRENTS 

Equipped with Eq. (H]), ([5]), and (E]), we are ready to investigate the surface bound states 
and spin currents at the surface. We do so by evaluating the spin and momentum- resolved 
densities of states along different spin projection directions. For example, for a given mo- 
mentum, the density of states at the surface for spins along positive (negative) z axis is 
given by p^±{e) = -^ImTr[i^g^(k, e; x = 0)], where g^ = g|i,,^^e+i5,£„>o is the retarded 
Green's function. Similar formulas apply by replacing z with the other corresponding di- 
rections. Here g is the 2x2 matrix in spin space given by the upper left block of g in the 
Nambu space. 

We shall show the results for representative cases of nontrivial Z2 with lA^/A^I = 2 and 
trivial Z2 with |Ap/As| = 1/2. Figj2] shows the spin-resolved density of states pz,± with 
0=±| and |Ap/As| = 2. The momentum ky = fcirsinc/) is a good quantum number as a 
result of the translational invariance along y. The peaks below the gap indicate the surface 
bound state energy corresponding to 



In the case of pure p-wave (A^ = 0), the bound state spectrum is given by Es^=±{ky) = 
=F|Ap| sin0.— i2i In this limit these surface quasiparticles have a fixed spin orientation. Hence 
Pz±{^) consists of a single delta function peak at e = =F|Ap| sin0, whereas Px,y±{^) are both 
identically zero. In contrast, in Figj2] (a) and (b), we find two subgap peaks, though with 
large differences in height. These indicate that the quasiparticles in the NcSC are no longer 
eigenstates of in the presence of order parameter mixing. These minor peaks in P2±(e) 
vanish eventually as |Ap|/|As| is increased toward infinity. Moreover, Px,?/±(e) (Fig [3l H]) 
are also finite at the bound state energies. This indicates that the quantization axes for 
the spins are not aligned with z. Furthermore, the spin resolved density of states are finite 
also for energies larger than |A_|. This contribution is a continuous (not delta) function 
in energy, indicating that this contribution to the spin densities and spin currents arises 




det {A(k),B(k)} 



0. 



(7) 
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from continuum states not bound near the surface. We note that kinks also appear in these 
density of states at the gap edges where e = ±|A-|-| = ±|Ap ± As|. In this case, they are 
±|As| and ±3|As|, respectively. 

For a given angle which parameterizes the quasiparticle path, the surface bounds states 
appear in pairs of equal but opposite energies. The bound state energy for the E > branch 
versus the angle is plotted in FigJS] It can be seen that the Andreev bound states are 
pushed toward the band edge by the s-wave pairing order parameter. We study in addition 
(FigJSl right panel) the spin polarization Si = P'^+'f"-- of the Andreev bound state versus 
the angle 0. The sum in the denominator is in fact independent of the direction i and equals 
the density of states. The S'j's in some sense describe the spin direction of the Andreev 
bounds state. With increasing Ap/Ag, Sz approaches a step function while approaches 
zero. Thus we can conclude that the effect of s-wave order parameter is to tilt the spin of 
Andreev bound states toward the x axis. 

The spin current density Jy at the surface is obtained via the expression, 

= 0) = ^Npvp J ^(sin0)T^Tr e„; 0)] , (8) 

The summation over the Matusbara frequency (T J2n Tr[..g]) can be replaced by the integral 

^ImTr[..g^(e)]tanh(^) 

with respect to the real frequency e. We shall first note an elementary symmetry relation 
followed from TRS of our superconducting state 



which is valid for all spin projections n. This can be seen in Figj2l Figj3]etc. From Fig El we 
see that we also have, in addition, the symmetry Pz±{ky, e) = Pz±{—ky, — e) (and consequently 
Pz+{ky,e) = Pz-{ky, —e)). The asymmetry between pz+{ky,e) and Pz-{ky,e) (and between 
Pz+{ky,e) and pz+{—ky,e)) causes a finite spin current Jy in the topological nontrivial cases 
as well as in the trivial cases, whereas eq guarantees that all spin accumulations are zero. 
For spin projections along x and y as in FigJ^land FigJU we have in addition p^y+{ky,e) = 
Px,y+(^?/, — e), (and similarly for H — > —) which forbids both the spin currents Jy and J|. 

In the trivial Z2 case where we choose |Ap| = |As|/2, no subgap state is found as in 
Figint^a). However, there is still an asymmetry between pz+{ky,e) and Pz-{ky,e) for the 
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continuum states between the two gaps, so that these states continue to contribute to Jy. 
The quasiparticles upon reflection at the surface is now analogous to transmission through a 
Josephson junction between two unequal-gap SC's with a relative phase difference,— where 
inter-gap continuum states contribute to a finite current. We suggest here that the contri- 
bution to the spin current can be pictured in similar manner. With Eq. (|H]), we find that Jy 
is indeed nonvanishing. This is justified in the real frequency domain as shown in Figj6]^b) 
where a finite contribution to Jy is found at energies e between the gaps — |A_|_| < e < — |A_| 
(|A_| < e < |A+|). Note that there is no contribution from states with |e| > |A+| (These 
two statements also hold for the topologically non-trivial case). 

In Figl?! the values of Jy{x = 0) as a function of the triplet to singlet order parameter 
ratio |Ap|/|As| is presented. Jy{x = 0) increases from zero in pure s-wave case toward the 
value obtained in the pure triplet case. In the topologically trivial regime, |Aj,| < \As\, 
Jy{x = 0) seems to be quadric in |Ap|; while in the topologically nontrivial regime, it is 
roughly linearly with |Ap|. The asymptotic value for Jy{x = 0) coincides with our previous 
results.— FigE] shows the values of the total spin current ly = dxJy{x). 

IV. DISCUSSIONS AND CONCLUSIONS 

The pairing Hamiltonian can be conveniently written in the helicity basis, H/^ = 
\ Sk(^A+e~*'^''a[^aLk+ ~ i/^-e'^'^^a\_a^_^_) + h.c. Here A-t are both real numbers so that 
the above is time-reversal invariant. 0k is the angle between the momentum k and the x 
axis. The relation A^p = (A+ ± A_)/2 can be obtained upon transforming back to the 
normal spin basis {ak|,ak^}. Thus the pure p-wave Rashba triplet SC and the pure s-wave 
singlet SC are recovered when A+ = =fA_., respectively; while the case with |A+| ^ |A_| 
belongs to the NcSC The nontrivial Z2 class corresponds to A+A_ < 0.— 1^ 

The pure s and p-wave cases, |A+| = |A_|, are well understood from the previous 
papers.—"^ The main difference between the pure triplet case and the present NcSC is that 
the surface Andreev bound states do not have a fixed spin projection. Moreover, besides 
the surface bound states, the continuum states between the two gaps, |A-|-|, also contribute 
to the spin current and are solely responsible for its nonvanishing value in the topologically 
trivial case. 

In our calculations, where the normal state spin splitting is not included, only the com- 
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ponent Jy is nonvanishing while the other components do not appear, even though the 
quasiparticles have their spins not ahgned with the z axis. Mixed-parity SC order param- 
eter along with the presence of Rashba spin-orbital interaction accounting for the absence 
of inversion symmetry may be closer to the physically accessible regime/^^ However, the 
inclusion can in principle generate nonzero spin currents even in the normal state under 
equilibrium.— Even so, the spin current obtained in the normal state for a spin-split band 
is much smaller than what we have obtained in our mixed-parity superconducting state 
without the corresponding spin-orbital coupling, provided that the factor TcEp/a^ is suffi- 
ciently large,— which is presumably true in usual situation when the spin splitting energy 
a <^ Ep. Moreover, inclusion of the spin-orbital coupling in the normal state will make 
the spin current ill-defined because of the lack of spin conservation. In Ref.— , the authors 
define the real spin current in superconductors by subtracting the contribution obtained in 
the normal state. They then obtained also finite Jy and J^. 

In conclusion, we use the exploding-decay tricks to obtain the quasi-classical Green's 
functions associated with a singlet-triplet mixed noncentrosymmetric superconductor. For 
the topologically nontrivial NcSC, we obtain a pair of Andreev bound states without a fixed 
spin projection and a consequent spin current Jy. For the topologically trivial NcSC, a finite 
spin current Jy remains even though no Andreev bound state can be found. 
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Appendix A: exploding & decaying trick in NcSC 

Here we shall present our scheme for obtaining the solutions given in Eq. (jlj) to the 
Eilenberger equations Eq. ([T]) and ([2]) for the NcSC in the presence of a boundary as shown 
in Fig{Tl Eq. (j4]) is expressed in terms of the exploding and decaying solutions Eq. ([5]) and 
(j6]) of an auxiliary problem that is identical to the present NcSC for x > 0. To make our 
reasoning clearer, we first consider the special case where d coincides with z, then eventually 
obtain the Green's function for general d. (A similar method to obtain the quasiclassical 
Green's function for a uniform non-centrosymmetric superconductor has been also used in 
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Ref.— . Our treatment extends to the non-uniform case and we also point out some useful 
mathematical relations not noted in there). The matrix M = (ze„ — A) appearing in the 
commutator in Eq. ([T]) becomes, after an obvious rearrangement of rows and columns, block- 
diagonalized. Explicitly, with m = k ■ R, it has the form 




+ tVFdu I 1=0, (Al) 



where the diagonal elements of M are 



M^^= ' " I • (A2) 




The parameters A± = ± Ap. By imposing the spatial dependence of e to g, Eq. 
(lAip becomes a problem for finding eigenvalues A and the corresponding eigenvectors. For 
example the upper diagonal block (++) simplifies to [M^^,g^^] ~ 2iXg'^'^ = 0. This re- 
lation is analogous to the pure s or p wave case and can be solved in the same manner. 



The eigenvalues A form the set {0,0, q;+, — «+} with a+ = ^ e\ + jA+P and the associated 
eigenvectors are respectively the set of 2x2 matrices {— ivrl, (7++, a^+, 6]*]^^^^}. The previous 
two are the "constant solutions" while a^^ and 61^^ are the decaying and exploding ones, 
given by 



-"^^'^ ,-^t'"rr"'i. (A3) 



-Alfa, -e„) ilA, 



12 



i|A+|2 -A+(a+-e„) 



12 



-A;(a+ + e„) -i|A 
respectively. Besides, the following equality holds as well as in e.g. Refj 



(A4) 



23 



= -i7r[a++,6++]{a++,6++}-^ (A5) 



— vr. 



(A6) 

«+ 

where the subscripts for eigenvalues are omitted if no confusion is generated. Similar results 



can be obtained in the ( ) block by defining «_ = ye^ _|_ |A_p. As for the off-diagonal 
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blocks (H — ) and ( — h), they share the same set of eigenvalues, {±as,±«d}, with the sub- 
scripts standing for the sum and difference by = . The eigenvectors belonging to 
the block (ij) = (H — ) or ( — h) will be labelled as {a^^^,^^^^}- It should be reminded that 
the decaying (exploding) solutions are associated with positive (negative) eigenvalues. 

Since the product of two solutions to Eq. (lAip also solves that equation, this product 
must be proportional to another solution, or it must vanish, depending on whether the sum 
of the corresponding eigenvalues coincides or not with any of the allowed eigenvalues. More 
precisely, if g^^ is a solution in the (ij) block with eigenvalue A, then the product is 
identical to Q^^^j^y^^ (up to a proportionality constant) or zero depending on whether Ai + A2 
coincides with one of the eigenvalues associated with the block {ik). Using Eq. (1A5I) . the 
following relations can be easily shown. 



(A7) 
(A8) 
(A9) 
(AlO) 
(All) 
(A12) 

Next we return to general d direction. Observing that g^"*^ in Eq. flASP can in fact be 
written as [— ^ii|^M((i = z)] in the 4x4 Nambu notation, we see that the information 
about d is embedded in the projection operator P± = l±^£i2.^ where a = {(Jx, cTyT^, a^},^ 
and in the order parameter A which appears inside M. From this, we obtain eqs (jHED- 
Indeed, it is useful to note that P± commute with M = (ie„ — — Ap). It is easy to 
see that P-|-M are the homogeneous solutions to eq ([1]). We can also see that the decaying 
and exploding solutions in Eq. ([5]) and ([6]) do satisfy eq ([1]), and one can further convince 
himself /herself by checking the relations, (a++)^ = (6++)^ = (a )^ = (b )^ = 0, and 
(gfe)^ = (g^"*" + g;^ )^ = "'""^j again hold. By similar reasoning as the d = z case, Eq. 
fR7l) - flA12p in the 4x4 form still hold for general d. 

Now we show how the present problem in FigH] can be solved in terms of these auxiliary 
solutions. For the outgoing path denoted by k, the most general solution with the correct 
limit as X — )■ 00 is 
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g{u) = g,(k) + c++(u)a++(k) + cr-(u)a-(k) + c+-(u)a+-(k) + cr+(u)a;+(k) , (A13) 

with u = K ■ k positive. For simplicity, the spatial dependence e~^"+"/*'^ associated with 
has been absorbed in the c- number coefficient c^^ , and similarly for the others. For 
the incoming path, 



giu) = a(k) + c++(u)b++(k) + C2 -(u)b-(k) + c+-(u)bi-(k) + C2 +(u)b-+(k) , (A14) 

with M = R ■ k negative here. We should remark that in the above equations decay- 
ing/exploding solutions associated with eigenvalue in the off-diagonal blocks has been 
excluded. We can see this by two different arguments. First, we can imagine that, far away 
from the surface, the inversion symmetry is restored. In that case — )■ and hence 
they would become constant solutions, which cannot exist at a; — )■ oo. Second, using 
relations eq( 1A7tlAT2|) . we can check that their appearance would violate the normaliza- 



tion condition ([2]). Now putting u = in Eq. flAlSP and multiplying this equation with 
i(k) = pia++ + paa— + p3a+- + P4a-+, and using Eq. (|A7l), (jM]), and (lATOj) . yields, 

i(k)g(0) = i7rA(k) , (A15) 

where the c-numbers {pi] i = 1,2, 3, 4} are arbitrary. Similarly, multiplying -B(k) = qib^'*' + 
qab"" + q3b+" + q4b^+ with Eq. f lATill and using Eq. (jM]), f lATTll . and f lAT2ll . yields 

E(k)g(0) = -i7rB(k) , (A16) 

where again {qfji = 1,2,3,4} are arbitrary. Multiplying the above two equations by -B(k) 
and A{k) respectively and add, one obtain the final expression in Eq. (jl]). The choices for pi 
and Qi are not restricted, as long as the anticommutator {A(k),B(k)} has a non- vanishing 
determinant and hence it is invertible. The determinant vanishes when the energy coincides 
with the Andreev bound state for a given 0, but the invertibility can still be guaranteed by 
adding a small imaginary number to the energy. Given two sets of matrix, Ai(k) and Bi(k), 
A2(k) and B2(k) will then yield the same g(0). Showing this is quite straightforward if one 
notices that ffist, the commutator and the anticommutator in Eq. (jl]) commute with each 
other, and second, Ai(k)A2(k) = Bi(k)B2(k) = 0. 

12 



To evaluate the spin current density at a general position x, we need the traces^^ 
Tr[o'r3^(?i)] at u. If d \\ z, we see that contributions to TilazT^g^u)] arises only from the 

++ and blocks, whereas TTla^T^glu)] and Tr[crj,^(M)] arise only from the H — and — h 

blocks. Hence the latter two have u dependence given by e-'^°'!>H/'"f (^^ee eq f lAlSp and 
(IA14p : note also Tilar-^gb] = 0). Hence for general d, if n is a vector perpendicular to d, 
Tr[(n • cr)T3g{u)] will have this same u dependence. Since for our superconductor d is in the 
x-y plane, Ti[a^r3g(u)] is simply Tr[cr^r3^(0)]e~^"'''"'/''^. The x integral needed for the total 
spin current is simply /q°° (ixe~^"''"'''^-f = ailp^^ 
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FIG. 1: The interface between vacuum(gray) and a noncentrosymmetric superconductor(white) 
with singlet-triplet mixed order parameter specified by text. The directions are defined by the 
shown axis. The quasicparticle is incident from the lower right of xy plane along the path denoted 
by k. (j) is the reflection angle between the x-axis and the outgoing path k. 




FIG. 2: (Color online) Momentum and spin-resolved surface densities of states /9z^± in unit of 
with |Ap|=2|As[. (a)(/) = ^. (b)(/) = Note that the upper(lower) plot in (a) is identical to the 
lower(upper) one in (b) due to the time-reversal symmetry. The numerical values associated with 
the subgap peaks are related to the small imaginary number 5 = 10~^ we used in the transformation 
ien — )• e + which is also true in FiglSl 
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FIG. 3: (Color online) Momentum and spin-resolved surface densities of states px,± in unit of ^ 
with |Ap|=2|As|. (a) (/) = ^. (b) (j) = Besides the same symmetry followed by TRS, it shows 
in additional symmetry with e o — e which forbids the spin current Jy. 
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FIG. 4: (Color online) Same as Fig [3] except here that py^± is shown. The results for (p = it^ are 
identical. 
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FIG. 5: (Color online) Left: Positive energy branch of the surface Andreev bound states, corre- 
sponding to the zeros of the anticommutators in Eq. d?]), as a function of the angle (p for two ratios 
of Ap/Ag. The upper(green) and lower(purple) are for the ratios of 1.5 and 3.5, respectively. E is 
measured in terms of A_. As the triplet component becomes larger, the bound state behaves more 
like in the pure triplet case where E = \Ap\ sine/) (the black dashed line). Right: Spin polarization 
of Andreev bound states (defined in text) versus the angle (j). Solid lines are for the ratio Ap/Ag 
of 1.5 and the dashed ones are for 3.5. The direction of spin polarization evolves in the xz plane 
as (/> varies. For smaller angles or larger Ap/A^, the bound states are more spin-polarized along z 
as in the pure triplet case. 




FIG. 6: Topologically trivial NcSC with |Ap|=0.5|As|. (a): Spin-resolved surface density of states 
Pz± associated with ky = kp/^. No subgap states appear, (b): Non-vanishing inter-gap direction- 
integrated contributions /q'^^^ d(/)(sin (?!))Tr{cr2r3[— ^Im g^(e, </>)]} to Jy in the unit of hvFNp\As\. 
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FIG. 7: Surface spin current density Jy (in unit of fivpN^lApl) versus the order parameter ratio 
|Ap|/|As| at a temperature T = |As|/100. In the large |Ap| hmit, Jy = ^hvpNplApl, while the 
spin current vanishes in the pure singlet case. As |Ap| increases from the topologically trivial to 
nontrivial regimes, Jy follows different dependence on |Ap| around the transition |Ap| = 



0.25 




FIG. 8: Total surface spin current density ly (in the unit of h'^Npyp) versus the order parameter 
ratio |Ap|/|As| at a temperature T = |As|/100. 
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